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Abstract. The Renyi entanglement entropy is calculated exactly for mode- 
partitioned isolated systems such as the two-mode squeezed state and the multi- 
mode Silbey-Harris polaron ansatz state. Effective thermodynamic descriptions of 
the correlated partitions are constructed to present quantum information theory 
concepts in the language of thermodynamics. Boltzmann weights are obtained from 
the entanglement spectrum by deriving the exact relationship between an effective 
temperature and the physical entanglement parameters. The partition function of 
the resulting effective thermal theory can be obtained directly from the single-copy 
entanglement. 


1. Introduction 

Entanglement is a central topic in quantum information theory. It is viewed as a 
resource for manipulating quantum information and accomplishing quantum computing 
tasks. In an isolated pure state pab = It/’) (t/’l, entanglement is manifested as quantum 
correlations between its partitions A and B. That is, the results of measurements done 
on a partition, say A, are statistically distributed according to the reduced density 
operator pA = trsPAB, where the partial trace tr^/-) is a quantum average done over 
a complete set of states for the complementary partition B. For these bipartite pure 
states, Pa and ps have identical spectra (eigenvalues) and hence entanglement measures 
can be calculated from either reduced density operator. When the partitions A and 
B are unentangled both pA and ps are projectors onto a single state vector. On the 
other hand, when A and B are maximally entangled pa and ps correspond to uniform 
distributions over their respective eigenvectors. That is, there is more entanglement in 
Pab when there is less certainty in the state of a given partition when the effects of 
its complement are averaged over. A similar relationship between the uncertainty in 
knowing the (reduced) state of a subsystem and the strength of its coupling with its 
complement or environment is well-established in statistical mechanics [1]. 
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The main pedagogical objective in this article is the presentation of several key 
concepts of quantum information theory using the language of statistical mechanics 
and thermodynamics. The construction of effective thermal models for a partition is a 
very useful concept that underlies more advanced ideas such as the emergence of area 
laws in the scaling of entanglement entropy [2-4] and the replica trick of quantum held 
theory [5, 6 ]. The systems under consideration here are also particularly illustrative: For 
example, uncertainty and huctuations in a mode-partitioned squeezed state are identical 
to that of a harmonic oscillator at thermal equilibrium. Some of the results given in this 
paper have been demonstrated in the past [7, 8 ], but here we emphasize the consonance 
of the quantum information and statistical mechanics perspectives. Additionally, we also 
present some new exact results such as the detailed calculation of the Renyi entanglement 
entropy in bosonic two-mode squeezed states and a multi-mode Silbey-Harris state [9], 
and the identihcation of the single-copy entanglement [ 10 ] as a free energy analogue 
[11] for these systems. For a broader overview of the research on entanglement in 
many-particle systems, the reader may refer to extensive technical reviews found in the 
literature [12,13] as well as others at a more introductory level [14]. 

The mode-entangled examples studied here involve coherent (displaced) and 
squeezed boson modes that are encountered in, for example, quantum optical systems 
[15,16], ultracold coherent matter systems [17,18], and nanomechanical systems [19]. 
We partition these composite states with respect to the mode indices A and B and 
analytically obtain the reduced density operator for a partition p^. This allows us to 
calculate entanglement measures (section 2 ) for these bipartite pure states, such as the 
Renyi entanglement entropy [20] and entanglement spectrum [21]. The entanglement 
spectrum is further used to construct effective thermal models that yield identical 
fluctuations in observables that are local to a given partition (mode-local). We hrst 
consider several cases of two-mode coherent and squeezed states (section 3) and then 
proceed to the case of the multi-mode Silbey-Harris state (section 4). 

2. Quantifying entanglement 

The reduced density operator pA is easily obtained when the composite state vector 
\'ip )ab normalized and expressed in Schmidt decomposed form [ 22 ]: 

d 

I^)ab = with {m\^\n)y = (1) 

n=\ 

and \n)^ is a Schmidt state vector with a subscript x G {A, R} that labels the mode 
partition, 7 ^ 0 is the corresponding Schmidt coefficient, and d is the Schmidt number 
or Schmidt rank of \'f)AB- Schmidt coefficients satisfy the normalization condition 
Yln=i\^n\‘^ = 1. For brevity, we will omit the mode indices for the state vectors 
\n)^ = \n) and the tensor product operator 0 when context is sufficient. Evaluating the 
partial trace over mode B then yields a diagonal representation for pA in the orthonormal 
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basis of mode A Schmidt vectors: 

d 

Pa = trs = '^Pn \n)A (^U > (2) 

n=l 

which is an operator of rank d. The eigenvalues of this reduced density operator 
are the complex squares of the Schmidt coefficients Pn = |A„|^. The state \'f)AB 
separable (unentangled) with respect to the bipartition AU B ii and only if there is 
exactly one non-zero Schmidt coefficient (or, equivalently, the Schmidt rank is d = 1). 
Furthermore, if there is more than one non-zero Schmidt coefficient (d > 1) then the 
state is entangled. For the special case in which |A„|^ = 1/d for all n G {1, 2,..., d}, the 
state is maximally entangled. Quantitative measures of entanglement in the bipartite 
pure state {'^\ab^ such as the entanglement spectrum and entanglement entropy, 

are discussed in more detail below. 

2.1. Entanglement spectrum 

The eigenvalues of the reduced density operator pA constitute a probability distribution 
over the Schmidt state projectors of the reduced system since = 1- This 

set of eigenvalues is referred to as the entanglement spectrum [21] and may be used to 
construct an effective thermodynamic (canonical) ensemble over the Schmidt vectors of 
A through 

Pa = --, Z = tiexpy-IJ ^En\n) {n\j. (3) 

n 

This identihcation allows one to dehne an effective reciprocal temperature fd and a 
so-called entanglement Hamiltonian H = 'Ylin^ri\n) {n\. We emphasize that this is 
an effective thermodynamic analogy, since the composite pure state has no dehned 
temperature (it is isolated). However, such effective approaches have been used to 
investigate the mechanisms behind entanglement generation, area laws, and some critical 
phenomena [21,23-27]. 

2.2. Entanglement entropy 

Commonly used measures for entanglement between partitions are the von Neumann 
and Renyi entanglement entropies, which quantify the degree of uncertainty associated 
with a density operator. When calculated for a reduced density operator, the Renyi 
entropy is a scalar function of the entanglement spectrum 

1 ^ 

S^^ = ( 4 ) 

1 fl 

' n=l 

where /x > 0 is the Renyi parameter. The von Neumann entropy Si is the special case 
Si = — YlnPn which is mathematically identical to the Gibbs entropy of statistical 
mechanics and the Shannon entropy of classical information theory. Other limiting 
values of the Renyi entropy include (1) the logarithm of the Schmidt rank = Ind, 
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(2) the Renyi-2 entropy S^^ 2 , which is related to the purity 7 = 'ffnPn — [28; 29], 

and (3) the single-copy entanglement (SCE), which quantihes the maximum 

amount of entanglement that can be distilled from one copy of a state using partition- 
local operations and classical communication [10,30-32], 


3. Coherent and sqneezed states 

We now turn our attention to the calculation of the preceding entanglement measures 
for several mode entangled coherent and squeezed states. Emphasis is made on the 
effective thermal theories for the reduced systems. 


3.1. Two-mode coherent state 


A two-mode coherent state \oi,(3Y is obtained by applying a two-mode displacement 
operator V{a,/3) on the two-mode vacuum state |0)^ lO)^: 

|a,/3>' = V(a,P) |0>., 10)^ = e”-'-”'-*'’*’-'’-''|0)^ |0>j,, (5) 

where the superscript “c” on \a,/3y denotes “coherent.” The displacement operator 
can be written as a product of mode-local operators by the Baker-Campbell-Hausdorff 
(BCH) formula, which immediately gives the one-term Schmidt decomposition 

|q, py- ^ |o)^ « |0)^ = |a>^ . (6) 


where a single-mode coherent state is 


a) 


n=0 



n). 


(7) 


A similar one-term result is also obtained when an arbitrary number state is displaced, 
yielding a coherent number state V{a,(3) |m)^ \n)^. Thus, a two-mode coherent state 
has zero mode entanglement entropy and, like the full composite system, each mode 
partition is in a pure state. This example illustrates an important property of bipartite 
entanglement: it is unchanged by the action of partition-local unitary operations [ 20 ]. 


3.2. Two-mode squeezed vacuum state 

A two-mode squeezed vacuum state |; 2 )® (where the superscript “s” means “squeezed”) 
is produced by applying the squeezing operator iS( 2 :) = on the vacuum state 

|0)x|0)b: 

°° pine\p.Y.\^n 

S(V |0)b |0)b = E |n)B I")b . (8) 

n=0 

where = re'^ is the squeezing parameter. This is a generalization of single-mode 
quadrature squeezing [33] to the case of two boson modes [34,35]. Comparison with 
the dehning formula (1) reveals that the last equality is the Schmidt decomposition for 



Effective thermodynamics of isolated entangled squeezed and coherent states 


5 



Figure 1: The Renyi entanglement entropy for a 
mode-partitioned two-mode squeezed vacuum state increases 
linearly with squeezing magnitude r for large r. The 
Renyi entropy is the same for displaced squeezed states and 
squeezed coherent states with identical squeezing parameters. 


a two-mode squeezed vacuum state jz)®. The entanglement spectrum {pn} therefore 
consists of the probabilities 


Pn = 


tanh 


2n , 


cosh^ r 


(9) 


which are independent of the squeezing angle 9. The Renyi entropy is readily evaluated 
because the entanglement spectrum forms a geometric sequence: 



ln(l — tanh^^ r) -|- 2 /r In cosh r 

/i- 1 


( 10 ) 


We hnd that the Renyi entropy increases monotonically with increasing squeezing 
magnitude r, where for large squeezing magnitudes Sf^{r —)■ oo) ~ 2 /ir/(p — 1 ), it 
increases linearly with r (hgure 1). Thus, if the squeezing magnitude is tunable in 
an experimental realization of a two-mode squeezed state, then one can, in principle, 
generate arbitrarily large amounts of mode entanglement entropy. 

Taking the limit p —>■ 1 leads to the known result [7], which we now identify as the 
von Neumann entropy: 


S\ (r) = cosh^ r In cosh^ r — sinh^ r In sinh^ r. 


( 11 ) 


This result is valuable since the theoretical calculation of the total von Neumann entropy 
in bipartite bosonic Gaussian states reduces to the evaluation of *S'i(r) for two-mode 
squeezed states [ 8 ]. Meanwhile, the Renyi-2 entropy is Sl{r) = In cosh 2r so that 
the purity is 7 = sech2r. This expression for the Renyi-2 entropy is useful because 
it is an experimentally accessible quantity that is important for the measurement of 
entanglement in Gaussian states via homodyne detections [29]. Moreover, the single¬ 
copy entanglement is Sl^{r) = 2 In coshr. Since the SGE is a measure of the maximum 
entanglement that can be distilled from a given copy of a state, this exact result is 
relevant to experiments that carry out such distillations on squeezed states [36]. 
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Let us now construct an effective thermal model for the reduced single-mode state. 
We equate the elements of the entanglement spectrum with a corresponding effective 
Boltzmann weight 

tanh^'^ r 
cosh^ r Xln 

We then dehne an effective reciprocal temperature ft through = tanhr and 

obtain a harmonic oscillator entanglement Hamiltonian 

OO 

H = ^ hivn \n) {n\ , (13) 

n=0 



with an energy shift so that the ground state energy i?o = 0 (we justify this shift in 
the following paragraph). In other words, the statistical properties of the reduced state 
is the same as that of a harmonic oscillator with frequency u in thermal equilibrium; 
The complementary mode B interacts with A as though it were in a thermal bath at 
an r-dependent reciprocal temperature /3 and vice versa. 

Since the ground state energy was chosen to be i?o = 0, the largest element of 
the entanglement spectrum in the effective thermodynamic picture is po = '^/Z. Thus, 
taking the limit p —)■ cx) to get the SCE from the Renyi entropy gives 

= 21ncoshr = In^—= InZ. (14) 


Hence, the evaluation of the single-copy entanglement allows one to calculate the 
partition function and free energy of the effective thermodynamic model. In the case of 
the two-mode squeezed vacuum state, the effective free energy F = —Soo/fd is 


F = — ^ In cosh r = ln(l — e ^^), 

fj [J 

which is identical to that of a harmonic oscillator at thermal equilibrium. 


(15) 


3.3. Two-mode displaced squeezed and squeezed coherent states 

We now investigate mode entanglement in states that are both displaced and squeezed. 
The two-mode displacement and squeezing operators generally do not commute and so 
we consider these cases separately. 

A two-mode displaced squeezed vacuum state is denoted as 

\a,/3,zy" = V{a,/3)[S{z) |0)^ lO)^]. (16) 

The two-mode displacement operator is a product of partition-local operators and has 
no effect on the entanglement entropy. The Renyi entropy of a displaced squeezed state 
is therefore the same as that of the original squeezed state prior to the application of 
the displacement operator: 
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Next we consider entanglement in a sqneezed coherent state where the order of 
displacement and sqneezing is interchanged from the preceding example: 

\z,a,/3r = S{z)[V{a,/3)\0)^\0)^]. (18) 

Using the BCH formula to interchange the order of squeezing and displacement operators 
gives 

S{z)V{a, /3) = V{a, zl3*)V{zaf /3)S(z). (19) 

Thus, a squeezed coherent state is also a displaced squeezed state with the same 
squeezing parameter and the Renyi mode entropy for \z,a,(3Y^ also reduces to that 
of a squeezed state: 

SY{z,a,Y) = SY{a,Y,z) = Sl{r). ( 20 ) 


4. Multi-mode Silbey—Harris state 


The multi-mode Silbey-Harris (SH) state is 


L 


It), ® If )b - It). 


-f)l 


( 21 ) 


which is a variational ansatz for a spin or qubit (q) coupled to an iV-mode harmonic 
oscillator or boson bath (B) via an interaction of the form + al) [9]- This 

interaction leads to entanglement that can be measured and used to locate possible 
quantum phase transitions [37]. 

The variational parameters of the SH state are the components fk of the vector f 
that give the displacements in an iV-mode coherent state: 


r N 


|f)B = exp 




'-k=l 


N 


|0)i|0)2---|0)^ = (^|/,) 


'k ■ 


( 22 ) 


k=l 


In this work, we choose the displacements fk to be real for simplicity. The single-mode 
coherent states \fk)k \—fk)k ^ot mutually orthogonal and so the SH state as 
written in (21) is not Schmidt decomposed with respect to the partitioning between the 
qubit and boson subspaces. However, calculation of the reduced density operator on the 
qubit space is straightforward with the knowledge of the matrix elements (real /) 


if\, l/>» = (23) 

if\i \-f)t = (24) 

and the fact that the trace of a tensor product of operators is equal to the product of 
the traces of each factor in the product. Evaluating the partial trace of |f^^)(f^^| over 
the complete set of boson number basis states yields a qubit density matrix in the |ti) 
basis 


Pq 


1 —c 

—c 1 


(25) 
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Figure 2: The Renyi entanglement entropy for a Silbey- 
Harris state (partitioned with respect to qubit and boson 
spaces) depends on all the displacement amplitudes through 
f • f. Maximum entanglement is achieved exponentially fast 
with increasing f • f. 


with the constant 


N 


C = 


n \-fk)k = exp(-2f • f). 


(26) 


k=l 


Therefore, the entanglement spectrum for the mode-partitioned Silbey-Harris state is 
± |} and the Renyi and von Neumann entropies are 




5P(f) = 


(f) = 


1 -p ■ 

l + c 


In 


(l + c)^ + (l-c)^ 


In 


l + c 


2 ^^ 

1 — c 


In 


1 — c 


(27) 

(28) 


2 2 2 2 ' 
respectively. We find that mode entanglement in the SH state is exponentially sensitive 
to the displacements fk (figure 2): If any of the fk is significantly different from zero 
the entanglement entropy quickly saturates to its maximum value In 2 . 

An effective entanglement Hamiltonian for the reduced qubit state can be 
constructed such that it has a zero lower energy and an energy gap A: 


// = A|-> (- 


(29) 


where |—) is the Schmidt vector with the (smaller) eigenvalue —f }• This entanglement 
Hamiltonian leads to the effective canonical partition function for a two-level system 

Z = l + e-^^. (30) 

Dehning the effective reciprocal temperature 

= tanhf • f, (31) 

equates each element of the entanglement spectrum with a corresponding Boltzmann 
weight of this effective thermal model. Again, correspondence between the single-copy 
entanglement and In Z for a two-level system with an energy gap A follows immediately 

2 


S™(f) = In 


_1 + exp(- 2 f-f) 


= ln(l + e-^^) = lnZ. 


(32) 



















9 


Effective thermodynamics of isolated entangled squeezed and coherent states 


The effective free energy F = —Soo//3 is 


F = - In 

/S 


■l + exp(-2 f-f)' 

-lln 

1 

2 


_1 


(33) 


5. Concluding Remarks 

We have derived exact expressions for the Renyi entanglement entropy in bosonic two¬ 
mode squeezed states and the Silbey-Harris polaron ansatz for a qubit interacting with 
several boson modes. These exact results allowed us to (1) measure the entanglement 
present due to the interactions between the mode partitions, (2) construct effective 
canonical distributions with the same probability distributions as the reduced density 
operators, (3) express the effective reciprocal temperatures in terms of the physical 
variables such as the squeezing magnitudes and displacement parameters, and (4) 
identify the single-copy entanglement as the logarithm of an effective partition function. 

In particular, we have demonstrated that in a two-mode squeezed vacuum state 
characterized by the squeezing parameter ; 2 , the Renyi entanglement entropy increases 
monotonically with r = \z\. Also, tracing out the effects of one boson mode leaves 
the other in a mixed state with identical statistical properties as that of a harmonic 
oscillator of frequency u at reciprocal temperature ft = — (Intanh^ r)/ha;. Increasing 
the squeezing magnitude increases entanglement and the effective temperature of the 
reduced subsystem. This mode entanglement is not affected by displacing either of the 
modes before or after squeezing. 

We further determined that entanglement between the qubit and a boson bath in 
a multi-mode Silbey-Harris state saturates exponentially to its maximum value of In 2 
when any of the boson modes are signihcantly displaced. In terms of the displacement 
parameter f, the reduced qubit subsystem is at an effective reciprocal temperature given 
by /d = — (Intanhf ■ f)/A, with A being the energy gap between the qubit states. 

Note added in proof: Effective thermal theories have also been constructed for 
entangled left- and right-moving bosonic modes in boundary states of some conformal 
held theories [38]. 
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